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I.  INTRODUCTION 

An  electromagnetic  emitter  is  defined  as  a source  which  emits 
electromagnetic  waves.  The  wave  as  it  propagates  out  from  the  source 
display  a variation  of  both  amplitude  and  phase.  The  points  on  the 
wave  which  have  the  same  phase  in  a given  region  may  be  thought  of 
joined  by  Imaginary  line.s,  known  as  isophasor  lines.  These  are 
analogous  to  contours  on  a map,  which  define  a constant  elevation.  For 
a point  source  emitting  electromagnetic  waves,  the  Isophasor  lines  would 
then  form  concentric  circles.  The  objective  of  this  report  is  to  design 
electromagnetic  emitters  whose  Isophasor  lines  would  not  form  concentric 
circles  with  the  center  of  tlie  circles  being  located  at  the  geometric 
center  of  the  electromagnetic  emitters.  This  implies  that  such  an  array 
of  electromagnetic  emitters  would  produce  over  a certain  angle  in  the  far 
field  an  electromagnetic  wave  whose  phase  center  may  not  coincide  with 
the  geometric  center  of  the  electromagnetic  emitters.  The  projected 
phase  center  may  even  be  outside  the  spatial  distribution  of  tlie  elec- 
tromagnetic emitters. 

The  report  is  organized  in  v.trlous  sections.  In  section  II  the 
design  of  skewed  isophasor  lines  is  presented  when  the  electromagnetic 
emitters  are  excited  by  a band-limited  signal.  The  analysis  has  been 
carried  out  in  the  frequency  domain.  Section  111  presents  the  anal- 
ysis of  skewed  isophasor  lines  when  the  emitters  are  excited  by 
arbitrary  time.  The  development  for  this  case  has  been  nresented 
in  the  time  domain.  In  botli  sections  II  and  III  the  relationship 
between  the  slope  of  the  isophasor  lines  to  tlie  measure  of  closeness 
between  the  projected  phase  center  and  the  geometric  center  of  the 
spatially  distributed  electromagnetic  emitters  has  been  developed. 

« 

»■ 
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In  section  IV  some  representative  computations  are  made  using  an 
optimization  method  due  to  Rosenbrock.  Appendix  A describes  the  salient 
features  of  the  optimization  procedure.  Appendix  B provides  the  computer 
program  listing  along  with  sample  input  and  output.  This  program  has 
been  utilized  in  obtaining  the  results  of  section  IV. 

There  are  various  applications  of  this  principle  of  design  of 
skewed  isophasor  lines.  One  such  application  is  to  the  phenomenon 
of  angle  noise  or  angle  glint.  Angle  noise  causes  a change  with 
time  in  the  apparent  location  of  the  target  with  respect  to  a refer- 
ence point  on  the  target.  The  apparent  angular  location  may  even 
fall  outside  the  target.  This  principle  may  thus  be  applied  in 
providing  misinformation  to  the  receiver  monitoring  the  target. 

This  is  because  target  locating  systems  can  be  divided  into  two  major 
categories:  phase  comparison  systems  and  amplitude  comparison  sys- 
tems. The  phase  comparison  technique  involves  two  or  more  receiving 
elements  separated  in  space  so  that  the  phase  of  the  received  signal 
at  the  two  or  more  points  can  be  compared.  The  phase  comparison 
systems  then  measures  the  tilt  of  the  phase  front  since  it  indicates 
the  target  to  be  in  a direction  normal  to  the  line  of  receiving 
elements  alignment  required  to  receive  the  signal  in  phase  at  the 
receiving  structure.  Thus  a phase  comparison  system  measures  the 
phase  front  of  the  received  signal  and  points  to  the  direction 
normal  to  the  isophasor  lines  of  the  received  signal. 

Amplitude  comparison  systems  generally  use  some  type  of  secondary 
receiving  elements  which  focuses  the  received  signal  to  a spot  in  the 
focal  plane.  All  amplitude  comparison  systems  by  some  means  or  other, 
locate  the  target  as  the  pointing  angle  of  the  receiving  elements  which 


center  the  spot  from  the  received  signal  In  the  focal  plane.  For 
example,  a conventional  amplitude  comparison  monopulse  radar  split  the 
spot  In  four  parts  by  a multiport  feed  such  that  the  spot  can  be  centered 
by  adjusting  the  receiving  elements  for  equal  amplitudes  In  each  port. 

The  significant  characteristic  of  the  amplitude  comparison  techniques 
Is  that  the  spot  location  Is  determined  by  the  Isophasor  lines  of  the 
received  signal.  Centering  the  spot  Is  accomplished  by  rotating  the 
receiving  structure  such  that  the  receiving  elements  align  with  the 
received  Isophasor  lines.  Consequently,  In  general  any  target  locating 
device  Is  essentially  a phase  front  measuring  device.  From  this 
measurement  the  device  tries  to  project  the  phase  center  of  the  received 
Isophasor  lines. 

Another  application  of  the  principle  of  designing  skewed  Isophasor 
lines  may  be  found  In  missile  guidance.  A missile  may  be  guided  by 
the  skewed  Isophasor  lines  and  thereby  reduce  the  complexity  of  the 
tracking  mechanism. 

A third  application  of  this  principle  Is  a very  Interesting  one. 
Goldman  has  shown  In  his  book  "Frequency  Analysis,  Modulation  and 
Noise",  that  by  the  application  of  the  principle  of  stationary  phase 
the  location  of  a signal  In  time  could  be  obtained.  Under  certain 
conditions  the  transmitted  signal  may  be  made  to  appear  originating 
from  a point  In  time  other  than  Its  true  location.  This  apparent 
origin  of  the  signal  may  even  be  outside  the  spatial  distribution 
of  the  electromagnetic  emitters,  as  the  time  Information  may  be  con- 
verted to  a range  Information. 

A fourth  possible  application  of  this  principle  may  be  in  Displaced 
phase  Center  Antenna  (DPCA) . It  may  be  possible  this  way  to  compensate 
for  the  deleterious  effects  of  the  airborne  radar  platform  motion. 
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II.  ANALYSIS  OF  SKEWED  ISOPHASOR  LINES  IN  FREQUENCY  DOMAIN 
2 . 1 For  CW  Waveforms 

When  the  transmitted  signal  from  a set  of  emitters  Is  CW, 
the  Isophasor  lines  of  such  a signal  is  easy  to  visualize. 

The  phase  front  In  this  case  Is  the  Isophasor  line  of  the 
set  of  emitters.  The  problem  then  is  to  select  a set  of 
emitters  which  would  produce  skewed  isophasor  lines  in  the 
far  field.  The  tilt  would  be  such  that  phase  front 
measuring  device  looking  at  the  waveform  would  project  a 
phase  center  which  may  or  may  not  coincide  with  the  geometric 
location  of  the  set  of  emitters.  In  this  section  a re- 
lationship is  also  derived  relating  the  tilt  of  the  skewed  j 

Isophasor  lines  to  the  distance  between  the  geometric 
center  and  the  projected  phase  center  of  the  spatially 
distributed  emitters. 

Consider  a set  of  emitters  situated  on  the  x - axis  ' 

as  shown  In  Fig.  1.  Each  emitter  is  excited  with  a complex  j 

1 

signal  amplitude  S (having  both  magnitude  and  phase)  and 
n 

they  are  separated  from  each  other  by  a distance  d.  The 
electric  far  field  pattern  due  to  2N  + 1 emitters  could 
be  obtained  as 

N I 

E = z S exp  ( sin  6)  (2.1) 

n A 

n— N 

where  X is  the  wavelength  of  transmission.  Let  the  complex 

excitation  amplitudes  S of  the  emitters  have  real  parts  a 1 

n n I 

and  imaginary  parts  as  b . Then  ! 
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^ ^ ^ ain  0)  + j sin  sin  e)|j 


..  (2.2) 


By  taking  the  sunnnatlon  now  for  only  positive  values  of  n 
In  eq.  (2.2)  leads  to 


+ {cos  (— sin  6)  + j sin  (— ^—  sin  0)} 


...  (2.3) 


Let  the  complex  far  field  E has  a magnitude  |e|  and  an  angle 
i(/.  Then 


. ,2rr  nd 
Bln  ( — — sin 


“ j Z {(a  + a ) cos  sin  6)  + (b  -b  ) 

[n-l  n -n  X -n  n 

- sin  0)}  4agj+  j ^bp  sin 


sin  0 


+ (b  +b  ) cos  (^■.-  *^  sin  6 ) } i 
n -n  A 


,..  (2.4) 


Hence  from  eq.  (2.4) 


tan  4' 


b + E Fca  -a  ) sin  sin  0)  + (b  +b  ) cos sin  0)1 

0 n-1  I " ^ ^ J 


r 

a )•  I I (a  +a  ) cos  ( 
^ n-1  t " 


2ir  sin  0)  + (b  -b  ) sin  ain  0) 

A -n  n A 


...  (2.5) 

The  phase  front  in  a particular  direction  0 has  a tilt 

d0.  By  carrying  out  the  derivative  operation  with  respect 
to  6 in  eq.  (2.5)  and  after  some  algebrlc  manipulations, 
the  tilt  is  obtained  as  follows: 


’1 


d^^  2iTd  Numerator 

"Tn  “ cos  0 r: 

d0  A Denominator 


(2.6) 


where. 


Numerator 


■f" 

n-1 


n{  (a  -a  ) cos  sin  0)  -(b  +b  ) 

n -n  A n -n 


sin  sin  0)} 


-(b^  ”*’-n)  sin  0)} 


a + Z { (a  +a 

] •[  i. 


\ / nd  j Q \ 

) cos  ( — — sin  0) 


n{  (a  +a  ) 
n -n 


sin  (^-T^"  sin  0)  + (b  -b  . cos  sin  0)} 

A n -n}  A 


sin  (-.—--  sin  0)  + (b  +b  ) cos  (— ;--  sin  0)}  + b 
A n —n  A 01 


][i. 
] 


{(a  -a  ) 
n — n 


...  (2.7) 


and. 

Denominator 


. /2Tr  nd  . 

) cos  ( — — sin  0) 


- (b^-b  ) sin  (^^^  sin 

n -n  A 


[ “O  “"n-l 

0)}  1 fb  + E {(a  -a 

J I 0 n-1  “ 
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sin  gjLn  0)  + (b  +b  ) cos  (IZLl!^  g^n  0)}  1 ^ 

A n — n A I 


...  (2.8) 


For  the  particular  case, 


N-1.  Sq  " 0.  bQ  - 0,  - 1,  a - a,  b^^  - 0.  b_^  - 0 


eq.  (2.6)  reduces  to 


^ ^ eos  0 

d0  X ” 


1 - a 


1 + a^  + 2a  (—  — sin  0) 


...(2.9) 


Incidentally  (2.9)  is  the  same  result  as  obtained  by  Meade  ^ij 

[4 


and  Dean  Howard 


Next,  given  the  tilt  ^ of  the  Isophasor  lines,  what 
is  the  separation  distance  between  the  geometric  center  of 
the  spatially  distributed  emitters  and  the  projected  phase 
center  from  the  phase  measuring  device.  A typical  set  up 
is  shown  in  Fig.  2 where  the  part  of  a skewed  isophasor 
line  PQ  is  shown  to  have  come  from  a spatially  distributed 
emitters  of  length  A.  The  isophasor  line  is  being  observed 
at  a distance  R which  is  in  the  far  field  of  the  emitters. 
Two  receiving  elements  P and  Q are  properly  aligned  along 
the  isophasor  line  so  that  their  measured  phase  difference 
is  zero.  So  the  projected  phase  center  would  be  along  TS 
which  is  perpendicular  to  PQ  and  is  tilted  at  an  angle  8 
to  OT.  The  object  is  to  find  the  length  L.  Now  the  path 
difference  between  OP  and  OQ  is  RQ(»^  dij)) . Here  dt|^  is  the 
phase  difference  in  the  electric  field  between  the  points 
P and  Q along  the  isophasor  line.  Also  PR  « R d0. 


Hence  tan  8 - ^ ^ ^ 

can  p 2ttr  ^0 
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• • • 


(2.10) 


Hence  L 


X d<l> 
2n  dd 


...(2.11) 


So  from  (2.6) 

. X dll)  . „ Numerator 

L ••  n d cos  0 j — — 

2tj  d0  Denominator 

Lm  Q Numerator  ...(2.12) 

A M Denominator 

Where  A • Total  spatial  width  of  the  emitter  distribution 

and  M - (Total  number  of  emitters  -1) 

Therefore  It  is  possible  to  make  the  ratio  ^ change 

with  the  complex  voltage  excitations  of  the  spatially 

distributed  emitters  and  their  Interelement  spacing  d . 

For  some  problems  It  Is  desirable  to  make  ^ arbitrarily  large 

A 

by  proper  choice  of  d , a and  b . In  section  IV  some 

n n 

examples  would  be  presented  to  Illustrate  these  points. 


2.2  For  Wideband  Waveforms; 

For  wideband  waveforms,  It  Is  very  difficult  If  not 
Impossible  to  define  Isophasor  lines.  But  It  Is  still 
possible  to  define  a projected  phase  center  for  a spatial 
distribution  of  emitters  transmitting  wideband  waveforms. 

In  order  to  obtain  such  a definition  for  a projected  phase 
center  It  Is  necessary  to  define  the  model  for  the  receiving 
elements  of  a wideband  waveform.  For  this  problem,  the 
phase  sensing  device  Is  assumed  of  the  form  shown  In  Fig.  3. 
Here  two  signals  x^^  *2  arriving  at 

the  two  phase  sensing  devices  A and  B at  a certain  instance 
t.  On  one  channel  the  Hilbert  Transform  of  the  signal  Is 
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performed  and  this  is  represented  by  of  the  original 
signal  and  then  the  two  signals  are  multiplied  and 
integrated  to  yield  the  output  I.  I is  then  related  to 
the  measure  of  closeness  of  the  projected  phase  center  to 
the  geometrical  center  of  the  set  of  emitters.  For  1=0, 
the  projected  phase  center  would  coincide  with  the  geometrical 
center  of  the  spatially  distributed  emitters  and  for  any  other 
values  of  I i*0,  the  two  centers  would  be  far  apart.  The 
distance  between  the  two  centers  could  be  obtained  by  taking 
known  waveforms  and  calibrating  I to  the  error  distance. 

For  the  special  case  of  a signal  originating  from  a 
single  source  Xj^  (t^  0^^)  and  (t,  ©2)  would  be  the  same. 

For  a real  signal,  the  signal  and  its  Hilbert  Transform 
are  always  orthogonal  and  hence  I would  be  zero. 
Mathematically, 


00 

I “ (t)  X (t)  dt  = 0 ...(2.13) 

00 

always  holds  if  x (t)  is  real.  So  for  a single  source 
transmitting  any  arbitrary  wideband  waveforms,  this  device 
would  make  no  error  at  all  in  prediction  of  the  phase  center 
of  the  spatially  distributed  emitters.  However,  for  a general 
case  of  multiple  sources,  the  situation  would  be  different. 
This  is  dealt  with  in  the  next  section  but  in  the  time 
domain. 
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III.  SKEWED  ISOPHASOR  DESIGN  FOR  ARBITRARY  TIME  SIGNALS 


For  this  case  a set  of  emitters  are  chosen  similar  to 
the  diagram  of  Fig.  1.  Instead  of  exciting  each  emitter  by 
a single  frequency  as  in  the  previous  case  of  2.1,  each 
emitter  is  now  excited  by  currents  of  the  form 

i (x,t)  “ A (x)  f (t)  ...(3.1) 

m m 

where  A (x)  is  the  spatial  distribution  of  the  current  on 
m 

emitter  m and  f (t)  is  the  time  dependence  of  the  current. 
For  an  array  of  point  sources  one  can  write 

N 

i (x,t)  = f (t)  Z A (x)  6(x-md)  ...(3.2) 

ni 


where  A (md)  represents  the  strength  of  the  discrete 
m 

emitter  at  x » md.  Now  the  far  field  at  a distance  R can 
be  expressed  as 


^ T^_  Z A F(t  - - + — tin  9) 


4ttR  „ m 
m“-N 


c c 


..(3.3) 


where  c is  the  velocity  of  light  and  is  the  permittivity 
of  free  space.  This  is  essentially  equation  (20)  of  Tseng 


Cheng  . 


The  model  of  the  receiver  is  assumed  to  be  same  as 


before.  This  is  Illustrated  in  Fig.  2.  So  for  this  case 


where  {•}  denotes  the  Hilbert  transform.  The  limits 
of  integration  correspond  to  the  times  when  the  Integrator 
starts  integrating  (T^)  and  when  the  Integrator  stops 
Integrating  (T2) . 

For  simplicity  it  is  assumed  that  f(t)  « sin  wt, 
so  that 

(t)}  ■ 0)  cos  (a)t  + 90°)  • -u)  sin  wt  ...(3.5) 

Hence , 


I 

t 

I 
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.^2 


J 8in{  2(u)t  - J)  + J (m  sin  9^  + k sin  62)  ^ dt 


2 , 


+ r 'sin{  -(m  sin  0,  - ^ sin  9„)}  ; dt 

j I ^ J 


> ...  (3.6) 


Since  the  purpose  of  the  Integrator  in  Fig.  2 is  a 
smoothing  operation,  the  average  value  of  the  first 
integral  would  be  zero  as  the  time  average  of  a "sin" 
function  is  zero. 

Let  02  * 6 

and  9^  - 9 -t-  A6  ...(3.7) 

where  A9  represent  the  azimuth  angle  subtended  by  the 

receiver  at  the  origin.  Then 
m sin  9j^  - k sin  62  - 


[V 


So 


t \ . 

5 ■ 2 ; -1  I sin  A9  , 

(m  cos  A0  - k)  + (m  sin  A0)  sln<^9  + tan  ^ ^ , 

...  (3.8) 

the  error  I can  be  expressed  as 


N N 


, . . ilil  i A, 


ojUo  2 ^ / \ 

I - -1/2  (^)2  E E A 1X2  - T 1 

m*-N  k—N  ^ I 

^ m 

Bin  ^ -^(m  cos  AB  - k)^  + gi^  sin  {e  + tan  | ° ^1-  _y| 

• « 

...(3.10) 

It  Is  Interesting  to  note  that  since  I Is  an  odd  function 
of  e. 


I (6)  d0  = A / sin  1 z sin 


(0  + ♦)  j 


r r 1 

J sin  U sin  (0  + i/;)J  d0  + j sln|i 


z sin  (0  + i|;)  I d0 


sin  z sin  (0  + III) 


,|  ae  - / 


sin  I z sin  (0  I d0 


...(3.11) 


This  Implies  that  the  total  error  Integrated  around  the 
set  of  emitters  Is  always  zero. 

Now  the  error  between  the  projected  phase  center 
and  the  geometrical  center  of  the  spatial  distribution 
of  apertures  would  be  given  by  (2.11). 


2n  d 


So  the  objective  now  Is  to  find  the  relationship  between 
I and  i|<.  Let 
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*1  ^*’1’®!^  ° 

and  “ sin  (cot  + tj^) 

then  from  a derivation  similar  to  that  presented  in  this 
section  it  can  be  shovm  that 


/ 


T -T 
2 1 

cos  cot  sin  (cot  + ip)  dt  « — x sin  ip 


. (3.12) 


Hence 

L „ -i-  dt{/ 

^ 2t7  de 


-If L. 

TT  t(T2  -T^ 


2 - 4i7 


d0 

...(3.13) 


So  it  is  possible  to  make  L larger  than  the  half  of  the 
spatial  distribution  of  the  emitters  by  proper  choice 
of  Che  amplitude  coefficient  of  Che  emitters. 

Some  examples  are  presented  in  the  next  section  to 
Illustrate  it. 
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IV.  EXAMPLES 


The  mathematical  formulations  presented  In  the  previous  sections  Is 
utilized  to  obtain  some  representative  numerical  values  for  various 
electromagnetic  emitter  configuration. 

Figure  4 represents  the  skewed  isophasor  lines  of  a 3 element  array. 
The  spacing  between  the  elements  is  1.5\  so  that  the  total  width  of  the 
aperture  Is  3. OX.  The  objective  then  is  to  obtain  the  largest  slope  for 
the  Isophasor  lines  within  an  angle  window  of  10  degrees  (between  175 
and  185  degrees) . The  variables  in  this  problem  are  the  complex  excita- 
tions of  the  elements.  The  phase  reference  for  this  diagram  is  referred 
to  the  geometric  center  of  the  array.  The  optimized  phase  of  the  far 
field  is  shown  by  the  solid  line  in  Figure  4.  The  dashed  line  represent 
the  gain  for  the  3 element  array.  The  gain  Is  defined  as  the  logarlthlmlc 
ratio  of  the  far  field  between  the  present  3 element  array  and  an 
omnidirectional  element  fed  with  the  same  power  as  the  3 element  array. 

It  Is  Interesting  to  note  that  within  the  angle  window  the  electric 
field  Intensity  remains  relatively  constant. 

Figure  5 represents  the  optimized  2 element  array  having  the  same 

aperture  length  of  3. OX  as  In  Figure  4.  This  represents  the  optimized 

case  of  Meade  [1]  and  Howard  [2].  It  Is  clear  from  this  diagram  that 

the  slope  of  the  phase  of  the  E-fleld  is  no  longer  linear  over  the  angle 

window  as  it  was  for  the  3 element  case.  The  gain  for  the  3 element 

case  was  higher  outside  the  angle  window.  An  Important  feature  Is  that 

the  linear  portion  of  the  slope  of  the  phase  Is  situated  near  a minimum 

field  Intensity.  Figure  6 presents  the  slopes  (2L)  for  the  2 and  3 element 

A 

cases. 
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Figures  7 and  8 represent  the  same  cases  as  presented  In  Figures 


4-6  except  now  the  aperture  length  Is  restricted  to  0.6X  Instead  of 
3 X-  Figure  7 represents  both  the  phase  of  the  far  field  and  the 
gain  of  the  3 element  array.  This  correspond  to  point  A on  Figure  10. 
Figure  8 represent  the  corresponding  quantities  for  the  2 element 
array.  For  Figures  7 and  8 there  Is  very  little  difference  both  in 
the  phase  of  the  far  field  and  the  gain.  Hence  the  criterion  (max- 
imize the  slope  of  the  far  field  by  proper  choice  of  excitations)  is 
not  a valid  one  for  closely  spaced  elements.  This  is  true  because 
for  closely  spaced  elements  the  array  becomes  very  sensitive  to 
excitations.  Hence  a more  meaningful  criterion  would  be  the  maximiza- 
tion of  the  slope  of  the  far  field  with  respect  to  the  excitations 
for  a fixed  source  norm.  This  implies  that  the  optimization  would  be 
carried  out  for  a fixed  power  input.  This  is  a very  interesting  area 
which  would  be  investigated  in  future  work. 

Figures  9,  10  and  11  present  results  obtained  for  various  windows 
and  minimum  interelement  spacing.  The  plots  are  the  ratio  of  the 
projected  phase  center  by  half  the  length  of  the  aperture  against  the 
number  of  elements  In  the  array.  The  maximum  slope  for  the  curves  over 
a given  window  has  been  obtained  by  varying  the  excitations  and  InCer- 
element  spaclngs.  Care  has  been  taken  In  the  optimization  procedure  to 
make  sure  that  the  Interelement  spaclngs  do  not  go  below  a certain 
value.  For  Figure  9,  the  minimum  Interclement  spacing  Is  0.5X«  It 
is  interesting  to  note  that  a large  ratio  ^ could  be  achieved  for  the 
same  number  of  elements  over  a smaller  window. 
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Figure  7 : Skewed  Isophasors  of  a 3 element  array 
having  aperture  length  of  0.6  X 
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Figure  10:  2L/A  for  various  elements  having  0.3  X as  an 
Interelement  spacing 
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Figures  10  and  11  display  the  same  characteristics  as  observed  In 


Figure  9 but  this  time  the  Interelement  spaclngs  are  greater  than  or 
equal  to  0.3X  and  0.2X  respectively. 

For  the  optimized  parameters.  It  Is  seen  from  Figures  9-11  that 
~ Is  always  greater  than  l.sX.  This  Implies  that  for  the  given 
windows  and  up  to  seven  number  of  elements  It  is  always  possible  to 
choose  excitations  which  would  produce  skewed  Isophasor  lines  In  the 
far  field.  These  Isophasor  lines  would  project  a phase  center  which 
Is  at  least  0.5  times  the  aperture  length  outside  the  physical  aperture. 

The  slopes  of  the  Isophasor  lines  In  the  far  field  becomes  much 
larger  as  the  elements  come  closer  together.  This  Is  seen  by  comparing 
the  curves  of  Figure  11  which  correspond  to  a d^  O.IA  to  those  of 
Figure  9 and  Figure  10  which  correspond  to  a d 2 0.3X  and  d ^ 0.5A 
respectively.  Hence  there  Is  a strong  Indication  that  It  may  be  possible 
to  achieve  larger  slopes  of  the  Isophasor  lines  from  a continuous  spatial 
distribution  of  electromagnetic  emitters  rather  than  a discrete  one. 

But  as  outlined  before  this  criterion  shall  be  used  along  with  an  addit- 
ional constraint  on  the  source  norm  for  the  arrays.  This  could  be  an 
worthwhile  area  to  look  Into  In  future  work. 
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V.  CONCLUSION 

The  possibility  of  obtaining  skewed  Isophasor  lines  from  two  emit- 
ters has  been  illustrated  by  Meade  [1]  and  Howard  [2].  From  their 
results  it  is  seen  that  large  slopes  of  the  Isophasor  lines  could  only 
be  achieved  over  a narrow  window.  The  theory  developed  in  this  report 
could  be  applied  to  obtain  skewed  isophasor  lines  over  large  angle 
windows.  Also  of  Importance  is  the  fact  that  sometimes  the  power  density 
is  considerably  greater  than  that  Indicated  by  Meade  [1]  and  Howard  [2]. 

The  theory  presented  in  this  report  can  also  be  applied  to  time 
waveforms  where  the  apparent  time  of  arrival  of  the  waveform  can  be 
changed  with  appropriate  phasing  of  the  electromagnetic  emitters. 

Another  Important  application  is  in  the  guidance  of  a missile  by  the 
skewed  isophasor  lines.  It  may  also  be  applied  in  the  case  of  Dis- 
placed Phase  Center  Antenna  (DPCA)  to  counteract  the  linear  motion  of 
the  airborne  radar.  However  the  technique  presented  here  can  find 
applications  to  many  other  electronic  systems. 
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APPENDIX  A - ROSENBROCK'S  OPTIMIZATION  METHOD 


The  function  - minimization  algorithm  due  to  Rosenbrock  [4^  la  well 
known  but  Is  Included  here  for  completeness. 

The  Rosenbrock  search  technique  uses  M mutually  orthogonal  directions 
during  each  search  cycle  to  find  a relative  minimum.  This  strategy  differs 
from  a steepest  descent  technique  which  uses  successive  orthogonal  directions, 
but  these  successive  directions  do  not  necessarily  form  a mutually  orthogonal  set. 
The  M mutually  orthogonal  directions  are  the  basis  for  the  success  of  this 
technique.  Moreover,  unlike  the  other  optimum  search  procedures  this 
method  does  not  require  any  derivatives  of  the  functions  to  be  minimized. 

The  basic  elements  of  the  Rosenbrock  algorithm  are  as  follows: 

a)  Step  Size:  The  step  size  In  a given  direction  Is  chosen  by  speci- 
fying an  arbitrary  magnitude  h and  then.  If  a step  h decreases  the  value  of 
the  function  (for  a minimization  problem),  h Is  multiplied  by  a constant 

(a>l).  If  the  value  of  the  function  Increases,  h Is  multiplied  by  a con- 
stant - 3(0<  3<  1).  h should  be  on  the  order  of  one  percent  of  the  average 
magnitude  of  the  variables.  To  determine  the  values  of  ot  and  3,  a series 
of  trials  using  various  values  of  a and  3 should  be  made  for  a given  class 
of  functions  for  which  the  function  to  be  minimized  belongs. 

b)  Direction:  The  Rosenbrock  algorithm  uses  M mutually  orthogonal  directions 
d^,  d2,  . . . , d^  at  each  stage  rather  than  choosing  a single  direction  In 

which  to  progress.  Therefore,  a search  Is  made  In  each  orthogonal  direction 
before  the  next  step  Is  chosen  (at  least,  one  trial  has  been  successful  - 
a value  less  than  or  equal  to  the  old  value  - and  one  has  failed  In  each 
direction) . 

There  are  three  cases  to  consider: 
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1)  The  first  trial  is  a success. 


ii)  The  first  trial  is  a failure  and  the  second  trial  is  a success, 
ill)  The  first  trial  is  a failure  and  the  second  trial  is  a failure. 

For  a start,  let  (identity  matrix)  where  d^,  d^,  ....  d^ 

are  individual  directions  in  which  steps  are  to  be  taken.  The  search  technique 

progresses  as  follows  (for  one  variable  x^^): 

1)  A step  is  taken  in  the  direction  d^  from  an  initial  point  x^ 

X » h (x  • d^) . 

2)  If  f (x^)  <f  (x®),  the  step  is  successful  and  the  step  size  is 
multiplied  by  0i  again  and  again  until  a failure  is  recorded.  After  this 
failure,  the  result  of  the  last  succcessful  trial  and  the  number  of  success- 
ful trials  is  stored. 

3)  If  f (x^)  >f  (x^)  on  the  first  trial,  a failure  results  and  the 
sign  of  the  step  size  h is  changed  and  the  search  continues  as  in  (2) . 

4)  If  f (x^)  > f (x^)  and  f (x^)>  f (x^) , the  second  and  further  trials 
are  determined  by  multiplying  h by  minus  B until  a success  is  recorded. 

After  a search  has  been  made  for  each  variable  in  all  directions  d^,  a new 
[^D^3  determined.  Let  g^  be  the  algebrlc  sum  of  all  the  successful  steps 
hj^  in  the  direction  d^  . Then  the  first  element  d^  of  Is  chosen  as  the 

vector  sum  of  all  the  vectors  I 8^^?.  The  other  elements  of  are 

determined  using  the  Gram-Schmidt  orthogonallzatlon  procedure.  Let 

Ai  - ^2^2*"  ’ ■ ’ 


^2  “ 


g2d5+  . . . .+g„d2 


where  A^^  is  the  vector  Joining  the  initial  and  final  points  obtained  by  the 
sum  of  vectors  d^,  and  A2  is  the  sum  of  the  successes  achieved  in  the  directions 
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other  than  the  first  and  so  on.  Then  from  the  Gram-Schmldt  orthogonallzatlon 
procedure 


1 ®1 
1 IB 


l®ll 


®2  * ^2  " ^2  **1 


!i- 

^ ■‘l«2l 


where  T is  the  transpose 


®M  “ - 


M-1 

Z 

IB=1 


1 

m 


1 "m 

So  for  the  second  iteration  the  searches  are  made  in  M orthogonal  directions 
of  process  is  continued  until  a convergence  criterion  is  satisfied. 

The  computer  program  listing  for  this  procedure  is  listed  below. 


1070  SUBROUTINE  ROTATE<X»Nf VFfUf Hf VV»NFIRST) 

1080  dimension  d(39)  rwOO)  »u(30f  30)  »a(30»30)  »b(30)  »un(30»30)  »ns(40) 
1090  bb=0.5 
1100  33=3«0 

1110  ifCnfirst-l)  15»15»17 

1120  15  do  24  i=l»n 

1130  d(i)=0.0 

1140  do  25  J=l»n 

1150  25  u(ifj)=0,0 

1160  24  u(i»i)=1.0 

1170  VV=Vf(H) 

1180  17  do  1 J=lrn 
1190  ns(j)sO 


k' 
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1200  1=1 

f 1210  e=h 

1220  vvt=vv 

1230  2 do  3 i=l»n 

1240  3 b < i ) =x ( i ) +e#u ( J » i ) 

1250  dv=vf<b)-vv 

1260  vv=vv+dv 

1270  do  to  (4.5»26»27) rl 

1280  4 1=3 

1290  if(dv)  28»28»29 

1300  29  1=2 

1310  vv=vv-dv 

1320  e=-h 

1330  do  to  2 

1340  5 1=3 

1350  if(dv)  28f28f35 
1360  35  1=4 

1370  41  if(ns<J)+20)  44»44»47 

1380  47  e=-bb*e 

1390  ris< J)=ns< J)-l 

1400  do  to  2 

1410  28  0=a3«e 

1420  ris( J)=ris( J)+l 

1430  do  to  2 

1440  26  if(dv)  28,28»42 

1450  42  e=e/3a 

1460  vv=vv~dv 

1470  do  to  44 

1480  27  dv=vvt-vv 

1490  if (dv)41»41r44 

1500  44  do  45  i=lfn 

1510  45  x(i)=x(i)+e*u(J»i) 

1520  1 d(J)=e 

1530  do  9 m=l»n 

1540  c=0.0 

1550  do  6 i=l»n 

1560  2Z=0.0 

1570  do  7 J=ifi»ri 

1580  7 z2=zz+d( J) #u( J» i ) 

1590  if  (m-l)  32»32»8 
1600  8 mro=ffi-l 
1610  do  11  J=lfn)m 
1620  z=0. 

1630  do  12  k=lrn 
1640  v=0. 

, 1650  do  13  l=mfri 

1660  13  v=v+d(l)*u<lrk> 

1670  12  z=z+v»uri< Jfk) 

1680  a<irifJ>=z 
1690  11  continue 

BEST  AVAIIABIE  COBf 
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1700  do  14  J=lfm» 

1710  14  zz=zz-3(iii»J>#uri<  Jri  ) 

1720  32  b<i)=sz 
1730  6 c=c+zziltzz 
1740  c=l»/sart<c) 

1750  do  10  i=lfn 
1760  10  un<m»i)=c#b(i) 

1770  9 continue 
1780  do  18  i=lfn 
1790  b(i)=0.0 
1800  do  19  J=l»n 
1810  19  u( i f J>=un( i » J) 

1820  18  continue 
1830  return 
1840  end 


I 

I 
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APPENDIX  B - COMPUTER  PROGRAM  DESCRIPTION 

The  user-oriented  computer  program  consists  of  a main  program  which 
supplies  all  the  necessary  data  required  for  the  execution  of  the  computer 
program  and  subprograms.  The  first  subprogram  FUN  In  this  case,  defines 
the  function  to  be  minimized  and  the  second  subprogram  ROTATE  is  the 
optimization  method  due  to  Rosenbrock.  The  latter  Is  described  In 
Appendix  A. 

The  main  program  supplies  all  the  Input  data.  The  Input  data  con- 
sist of  the  following  statements.  Statement  number  40  provides  the  value 
for  N (as  shown  In  Fig.  1)  through  MNE. 

Statement  numbers  90-120  provide  the  Initial  guess  for  the  excita- 
tions of  the  electromagnetic  emitters  In  the  first  N8  elements  of  the 
array  Y.  The  elements  of  Y(I)  are  defined  In  the  following  manner: 

- 0.5  [y(1)  + Y(2)] 

b^  - 0.5  [y(3)  + Y(4)] 

a_^  - 0.5  [y(2)  + Y(l)] 

b.i  - 0.5  [y(3)  - Y(4)] 

and  so  on.  The  statement  number  130  provides  the  Initial  guess  for  the 

Inter-emltter  spacing  for  the  array  through  Y(N9). 

2L 

Statement  number  150  defines  the  value  of  ^ through  the  constant 

SLOPE. 

Statement  numbers  170-180  define  the  angle  iflndow  over  which  the 
desired  slope  Is  to  be  achieved.  MP  provides  the  start  of  the  angle 
window  In  degrees  and  MQ  Is  the  last  angle  of  the  window,  taken  at  steps 


40 


of  MR  degrees.  It  is  not  necessary  to  have  the  window  synmetrlcally 
placed  around  the  z>axls  of  Fig.  1. 


Statement  number  250  specifies  the  maximum  number  of  Iterations  to 
be  made  by  the  optimization  method  to  yield  the  final  result. 

Statements  390,  400  and  480  provide  the  excitations  for  the  array 
elements  and  490  prints  the  value  for  the  Interelement  spacing  and  the 
value  of  the  function  through  VFVAL  at  the  end  of  each  Iteration. 

580  prints  the  slope,  the  magnitude  and  phase  for  the  far  field 
within  the  specified  window.  This  Is  done  at  the  end  of  the  final  Ite- 
ration. 

Statement  650  defines  the  constraint  Imposed  on  the  Inter-element 
spacing.  For  this  problem  Interelement  spacing  cannot  be  less  than  or 
equal  to  0.4  X . 

The  program  as  presented  represents  a 3 element  array  for  MNE-1. 

So  for  a 2-element  array  the  necessary  modifications  of  the  program  are 
as  follows.  The  statements  450-490,  870-880  should  be  deleted.  State- 
stent  80  should  read  N9  « N7.  Statement  641  should  be  Introduced  to  read 
N9  ■ N7.  In  statement  650,  the  contralnt  should  be  changed  to  half  the 
Interelement  spacing  desired.  For  example.  If  the  problem  be  such  that 
the  Interelement  spacing  for  a 2 element  array  be  restricted  to  be  greater 
than  or  equal  to  0.5  A then  Y(N9)  ^0.25 

The  computer  program  listing  is  presented  below. 


10  EXTERNAL  FUN 

20  COHNON  /AJ/MPfMQfMR»NS»SL0PE»MNE»N7fN8fN9»CC 
30  DIMENSION  Y(30) f AXES(30>30) 

40  NNE»3 
SO  NA>4*MNE 
60  N7>N61-1 
70  N8>N74^1 


80  N9-N8-I-1 
90  Y(1)»1*0  • 
100  Y(2)-0.01 


esi  m\m  copy 
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DO  11  I=3»N8 
11  Y(I)=0. 

Y(N9)=0.7 

H=0.01 

SL0FE=16 

PF=180./3, 14159263 

MP=177 

MQ=183 

MR=1 

MS=0 

J=1 

NQ=N9 

PRINT  189»MNE»MP»HQ»SL0PE»Y(N9) 

189  FORMATCO  ■ »3I9»2F10.5) 

DO  2 11=1 » 10 

CALL  ROTATE  ( Y»NQ»FUNf AXES»H» VFVALr J) 

NN=1 

CC=0 

DO  23  NL=1»MNE 
A1=0,5*(Y(NN)+Y(NN+1) ) 
Bl=0.5*(Y(NN+2)+Y(NN+3) ) 
A2=0.5»(Y(NN+1)-Y(NN> ) 
B2=0.5»<Y(NN+2)-Y<NN+3) ) 
C1=SQRT<A1»A1+B1*B1 ) 
C2=SQRT(A2*A2+B2*B2) 
D1=PP*ATAN2(B1»A1) 

D2=PP*ATAN2<B2»A2) 

CC=CC+C1*C1+C2*C2 
PRINT  108rNLrAlfBl»ClfDl 
PRINT  108»-NLf A2»B2»C2rD2 
108  FORMATC  •rI5»4F15.5) 

NN=NN+4 
23  CONTINUE 
NL=0 

C1=SQRT(Y<N7)*Y(N7>+Y(N8)»Y<N8)) 
D1=PP*ATAN2(Y(N8) fY<N7) ) 
CC=SQRT(CC+C1*C1) 

PRINT  108»NL»Y(N7)»Y<N8)fClfDl 
PRINT  100fY(N9) rVFVAL 
IF<VFVAL.EQ.O. ) GO  TO  5 
100  FORMATC  ■r4F15.7) 

J=J+1 

2 CONTINUE 

5 MS=1 
MP=MP-3 
Ma=MQ+3 
MR=1 

ZZ=FUN<Y) 

6 CONTINUE 
STOP 

END 


BBiAvwimtoe 
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620  FUNCTION  FUN(Y) 

630  niMENSION  Y(30> 

A40  COMMON  /A J/MP » MO » MR » MS » SLOPE » MNE » N7» N8 » N9» CC 
650  IF <Y(N9) .LE.0.4)  GO  TO  4 
660  DO  12  1=1. N9 

670  12  IP(ADS(Y<I)).GT.10.)G0  TO  4 
680  FLJN=0. 

690  DO  1 I=MP.Ma.MR 

700  PI=0.0174533*FL0Ar< I ) 

710  PH1=6.283186»SIN<PI )*Y(N9) 

720  UP=0.0 
730  DN=0.0 
740  NN=1 
750  U1=0. 

760  03=0. 

770  DO  2 N=1.MNE 
780  THT=N*PHI 
790  CT=COS(THT> 

800  ST=SIN(THT) 

810  UP=UP+Y(NN)*ST+Y(NN+2)*CT 
820  DN=DN+Y<NN+1 )*CT-Y(NN+3)*ST 
830  01 =01 +N* ( Y ( NN ) *CT-Y ( NN+2 ) *ST ) 

840  03=03+N*(Y(NN+l)»ST+Y<NN+3)»CT) 

850  NN=NN+4 
860  2 CONTINOE 
870  0P=0P+Y(N8) 

880  DN=DNfY(N7) 

890  SP=COS ( PI ) * ( 01 *DN+U3*0P ) / ( OP*OP+DN*DN ) 

900  SP=SP/FLOAT(MNE) 

910  IF  (MS.EQ.l)  GO  TO  11 
920  SP=ABS(SP) 

930  IF (SP.LT. SLOPE)  FON=FON+SLOPE~SP 
940  GO  TO  1 

950  11  AAG=SaRT(OP*OP+DN»DN) 

960  AAG=AAG/CC 

970  AAG=20  lALOGlOCAAG) 

980  SR=ATAN2(0P.DN) 

990  SI=180, 73, 14159263#SR 
1000  PRINT  101 . I .AAG.SI .SP 
1010  101  FORMATC  ■.I5.3F15.3) 

1020  1 CONTINOE 
1030  RETORN 
1040  4 F0N=100. 

1050  RETORN 
1060  END 


WflVAB'-t  CflP't 
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APPENDIX  C - COMPLEX  EXCITATIONS  OF  THE  OPTIMIZED  ARRAY 


The  complex  excitations  of  the  electromagnetic  emitters  that  yield  a 
prescribed  slope  of  the  phase  of  the  far-field  pattern  are  presented  in 
this  section. 

The  complex  excitations  required  to  yield  the  far-field  pattern 
presented  in  Figure  4 are  as  follows: 
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similarly  for  the 

field  pattern  for  the 

Figure  5,  the  excitations 

are  as  follows: 
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For  Figure  7,  the 
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are: 
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The  corresponding 

excitations 

for  Figure 
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Finally  the  complex  excitations  required  for  a 7 element  array  to 
yield  point  B in  Figure  9 are  as  follows: 
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